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We report 63 Cu(2) spin-lattice relaxation rate measurements of Yl^CusOg.gs in magnetic fields 
from 2.1 T to 27.3 T obtained from 17 0(2,3) nuclear magnetic resonance spin-spin relaxation. For 
T < 120 K, the spin-lattice rate increases with increasing magnetic field. We identify this magnetic 
field dependence with the change in the low-energy spectral weight originating from d-wave pairing 
fluctuation corrections to the density of states. 



PACS numbers: 74.25.Nf, 74.40. +k, 74.72.Bk 

Nuclear magnetic resonance has played an important 
role in elucidating the nature of high-T c superconductiv- 
ity In most metals the nuclear spin-lattice relax- 
ation rate divided by the temperature, 1/T\T, is a con- 
stant proportional to the square of the density of states at 
the Fermi level, Nf- In the normal state of many high-T c 
superconductors an increase in 1/T\T of planar Cu with 
decreasing temperature has been attributed to antiferro- 
magnetic spin fluctuations ||. At lower temperatures, 
in the superconducting state, the rate of planar Cu de- 
creases strongly with decreasing temperature as the gap 
in the quasiparticle spectrum develops. The crossover 
from normal to superconducting behavior occurs around 
100 K, substantially above the transition temperature of 
optimally doped YBCO. We have investigated this cross- 
over experimentally and theoretically. We show that the 
cross-over can be understood quantitatively in terms of 
pairing fluctuation corrections to the spin lattice relax- 
ation rate in optimally doped YBCO. 

Because of their large anisotropy and small coherence 
lengths, the onset of superconductivity in high-T c mate- 
rials is preceded by the effects of strong superconducting 
fluctuations on the normal-state properties, including the 
specific heat [EL diamagnetism || , nuclear spin- lattice re- 
laxation rate |^-|^], and Pauli susceptibility [|7|,|i"o||. Here 
we report on the field dependence of \jT\T of planar 
copper, 63 Cu(2), in optimally doped YBCO. We find that 
below 120 K the relaxation rate increases with increasing 
field with a typical field scale of 10 T. We quantitatively 
account for this behavior in terms of pairing fluctuations 
with d-wave symmetry 

Our aligned powder sample of 30% - 40% 17 0-enriched 
YBa 2 Cu 3 6 .95 has been investigated previously [jio] O. 
Our measurements cover the temperature range 70 K to 
160 K over a wide range of magnetic fields, from 2.1 T to 
27.3 T. The crystal c-axis was aligned with the direction 
of the applied magnetic field, the z-axis. Low-field mag- 
netization data show a sharp transition at T c (0) = 92.5 K. 



In order to study planar copper nuclear spin-lattice re- 
laxation we take advantage of its direct effect on the 
17 0(2,3) NMR spin-spin relaxation which we can accu- 
rately measure using a Hahn echo sequence: 7t/2-t-7t- 
acquire. Our typical 7r/2 pulse lengths were 1.5 /is, ex- 
cept at 2.1 T where pulse lengths were 2.5 /is, giving us a 
bandwidth > 100 kHz. After the 7r/2 pulse the process- 
ing nuclear spins dephase because of variations in the 
z-component of the magnetic field in the sample. The 
dephasing from static processes is recovered after the it 
pulse, leaving the echo intensity to be determined pre- 
dominantly by copper spin-lattice relaxation, as has been 
recently demonstrated The 17 0(2,3) (1/2 <-> -1/2) 

resonance has a low frequency tail owing to oxygen de- 
ficiency in a small portion of the sample |l2| . Its effect 
on our measurements can be eliminated by performing a 
non-linear least squares fit in the frequency domain for 
each echo, a method similar to that of Keren et al. |l6| . 
The success of this procedure was checked by relaxation 
measurements performed on the first low frequency satel- 
lite of the 0(2,3) [0. The oxygen resonance is much 
narrower than that of copper (by a factor of 6 at 8.4 T) 
and thus 17 NMR is more favorable for our experiments. 
This is particularly true for the high field experiments, 
H Q > 15 T, performed in a Bitter magnet at the Na- 
tional High Magnetic Field Laboratory in Tallahassee, 
Florida. The measurements for H a < 14.8 T were ob- 
tained with superconducting magnets. The reliability of 
this technique for measuring T\ of 63 Cu(2) was tested by 
comparison with direct measurements of T± performed 
on the same sample. 

We extract T x of 63 Cu(2) from 17 0(2,3) spin-spin relax- 
ation data following the proposal of Walstedt and Chcong 
]l7[ that the dominant mechanism for spin-echo decay of 
^O is the copper spin-lattice coupling. The z-component 
fluctuating fields from copper nuclear spin flips are trans- 
ferred to the oxygen nuclei by Cu-0 nuclear dipolar in- 
teractions. To account for this process Recchia et al. EH] 
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derived an expression for the 17 spin echo height, M(r), 
as a function of pulse spacing r, 
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We performed a nonlinear least squares fit of the data 
to Eq. (|l|) in the range 50 /is < r < 350 jtxs, with Ti of 
63 Cu(2) as a fitting parameter. The sum was performed 
over all Cu neighbors in a radius of 12 A; r, is the Cu-0 
distance; 0, is the angle between the applied field and the 
Cu-0 axis; is T\ of the i^ 1 copper nucleus; 7 = 3/2 
is the copper nuclear spin; k is an enhancement factor 
due to the Cu-0 indirect coupling which we determine 
to be 1.57; and Tib. is the Redfield contribution to the 
rate. An example of the fit is presented in the inset to 
Fig. [TJ at 19 T and 95 K, and is compared with the mea- 
sured relaxation profile. The fit to Eq. (Q) is sufficiently 
accurate that we can rely on its systematic behavior. We 
have also compared our 1/T\T data with direct measure- 
ments of 63 Cu(2) spin-lattice relaxation rates taken from 
earlier work fl pT|JT8| for several magnetic fields, as shown 
in Fig. g. The measurement at 7.4 T was performed on 
our sample [fblf . 

Our results for \/T{T are presented in Fig. [lj Above 
120 K there is no discernible field dependence within ex- 
perimental accuracy of ±2%. However, near the peak 
in 1/T\T we find that the rate increases with increas- 
ing magnetic field. At 95 K the rates at 2.1 T and 27.3 
T differ by 17%. The peak in l/T{T versus T shifts to- 
ward lower temperature as the field increases and the 
rate drops sharply in the superconducting state, consis- 
tent with reduction of T c by the field JTo| . We show below 
that pairing fluctuations are in quantitative agreement 
with this behavior, and that a purely magnetic mecha- 
nism with a spin pseudogap is difficult to reconcile with 
the field scale. 

In underdoped materials the temperature dependence 
of the Knight shift, K(T), and the peak in 1/TiT has 
been associated with the opening of a spin pseudogap 
|p9| in the spin excitation spectrum below a temperature 
T* > 100 K. The temperature scale T* was suggested 
to be a rough measure of the pseudogap, with a mag- 
netic field scale of H* = k B T*/fi Bl > 140 T. This ex- 
ceeds by far the field scale of ~ 10 T that we observe 
in 1/T\T in our optimally doped sample. The large field 
scale, H* 3> 10 T, for a spin pseudogap is consistent with 
recent neutron scattering measurements that show that 
the resonance peak of optimally doped YBCO remains 
almost unaffected in a field of 11.5 T |BC(|. 
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FIG. 1. Spin-lattice relaxation rate of 63 Cu(2) in YBCO 
as a function of temperature for the following fields: 27.3 
T (X), 22.8 T (O), 19 T (y), 14.8 T (^), 8.4 T (©), 5.9 
T (A), 3.2 T (□), 2.1 T (X). Dashed lines are guides to 
the eye. The solid line is a fit to (TiT)" 1 oc T X /(T + T x ), 
T x — 103 K H. Inset: Spin-spin relaxation of 17 NMR 
at 19 T, T = 95 K and a fit to Eq. <JTJ) . 

In high-T c materials superconducting fluctuations are 
expected to have a significant effect on \jT\ near T c . 
Diamagnetic fluctuations do not play a role in our mea- 
surements of T\ since they alter the magnetic field mainly 
along the axis parallel to the applied field; only transverse 
fields contribute to relaxation of the ^-component of the 
nuclear spin. The pairing fluctuation contributions to 
the rate result from fluctuation corrections to the den- 
sity of states (DOS) and from the Maki-Thompson (MT) 
corrections to the local dynamical susceptibility. The 
corresponding Feynman diagrams for these corrections 
are shown in Fig. ||. The propagators and vertices are 
defined below and in Ref. [Ql. A systematic analysis in 
T c jE F 1, where Ep is the Fermi energy, shows that 
all other pairing fluctuation contributions are negligible 
for local quantities like \jT\ ||. The pairing fluctuation 
correction is sensitive to the symmetry of the order pa- 
rameter fluctuations because of the difference in sign of 
the MT (positive) and DOS (negative) corrections, and 
because of the sensitivity of the non-s-wave pairing fluc- 
tuations to disorder. In the case of s-wave pairing fluctu- 
ations the dominant contributions to the rate come from 
the positive MT processes M, which are insensitive to 
non-magnetic disorder. A magnetic field suppresses the 
MT and DOS contributions, and leads to a suppression 
of the rate for s-wave. In the case of d-wave pairing the 
field dependence of 1 /T% is reversed compared to that for 
s-wave pairing. Scattering by non-magnetic disorder 
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FIG. 2. Pairing fluctuation corrections, to leading 
order in T c /Ef, for the nuclear spin-lattice relaxation 
rate, a) is the Maki- Thompson process, b) and c) are 
the density of states corrections to the rate. K is the 
impurity-renorrnalized pair fluctuation propagator. 

leads to strong suppression of the MT corrections for 
d-wave fluctuations. The DOS corrections survive non- 
magnetic scattering, but are suppressed by a magnetic 
field leading to an increase in 1/Ti with increasing field, 
even for modest levels of disorder. As we show below, our 
results provide a consistent and quantitative account of 
the field dependence of the nuclear spin-lattice relaxation 
rate above T c . 

In order to make a quantitative comparison between 
the leading order pairing fluctuation corrections and the 
experimental field dependence of the rate we isolate the 
fluctuation corrections to the experimental rate by writ- 
ing (TiT)tJj. = (TiT)- 1 + ^(TiT)- 1 , where the normal- 
state rate is fit to the AFM Fermi- liquid model ||, 
(TiT)" 1 cx T X /(T + T x ). We obtain T x = 103 K from 
a fit to high temperature data at 8.4 T. The fluctuation 
contributions are indicated by <5(TiT) _1 . These values, 
normalized by (TiT)" 1 , are plotted in Fig. || as a func- 
tion of magnetic field at 95 K along with our theoretical 
calculations of the pairing fluctuation corrections. 

The calculations of 1/Ti assume a quasi-2D cylindri- 
cal Fermi surface, with an isotropic in-plane Fermi ve- 
locity Vf. We expect the pairing fluctuations to be 
predominantly two-dimensional in a magnetic field be- 
cause of Landau-level quantization. A cross-over to 3D 
fluctuations is possible below ~ 2 T near T c . A sum- 
mary of the calculation is provided here; more details 
can be found in Ref. The pairing interaction is 

V(p, p') = r)(p) ■ g ■ r)(p'), where t](J)) is the normalized 
pairing amplitude; for s-wave pairing rj{p) — 1 while for 
d-wave pairing rj(p) = v / 2cos2?/>, where tp is the angle 
between the crystallographic d-axis and p. 

The pair fluctuation propagator is defined in terms 
of the sum over ladder diagrams in the particle-particle 
interaction channel; the propagator factorizes into 
VtpMQW), where T(Q) _1 = g' 1 -T^ B 2 (e n ,Q), 
B 2 (e n ,Q) = E P: V(P)V(P,Q)G(P)G(Q - P) and G(P) 
is the quasiparticlc Green's function. We use a short- 
hand notation: P = (e n ,p), P' = (e„/,p') for fermion 
quasiparticles, and Q = (u>i,q) for bosonic Matsubara 
energy and pair momentum of the fluctuation modes; 
the pair momentum, q, is quantized because of or- 
bital quantization in a magnetic field. We include dis- 
order via the standard averaging procedure for dilute 



impurity concentrations [ pT| . Impurity scattering in- 
troduces an elastic scattering time in the quasiparti- 
cle Green's function, G(P) = (ie n — S(e„) — £(p)) _1 , 
where = e(p) — p is the quasiparticle excitation 

energy, X(e„) = — (i/2r + i/2r^)sign(e„) is the self en- 
ergy, and r is the elastic scattering lifetime. We include 
inelastic scattering through the lifetime t^. Impurity 
scattering modifies the fluctuation propagator directly 
through a vertex correction in the particle-particle chan- 
nel, fj(P,Q) = r,(p) + E?v(p)G(P)G(Q - P)C(e„,Q), 
where C(e„,Q)- 1 = &- 1 - EpG(P)G{Q - P) is an 
impurity Cooperon-like propagator and a — 1/2tttNf 
is the impurity scattering vertex. The full impurity- 
renorrnalized pair propagator, K(P,P',Q), which enters 
the dynamical susceptibility diagrams shown in Fig. 0, is 
given by rj(P, Q)L{Q)i)(P' . Q). The leading order fluctu- 
ation correction to 1/Ti then follows from the Feynman 
rules for evaluating the diagrams and is given by, 

SxmM = -2\A | 2 B i(tn, Q)B 1 {e n - u m , Q)L(Q), 

n,Q 

S XD (u; m ) = 4|A | 2 J2 Gx{e n - ^n) ^ 2 ^ ) Q) T(Q), (2) 



(S(TiT)- 1 = lim 21m 

uj— >0 



(3) 



with Bi(e„,Q) = EpV(P,Q)G(P)G(Q - P), G x (e n ) = 

J2pG{P), and \A\ 2 are momentum-averaged hyperfine 
form factors M. We analytically continue Eqs. (0) to 



real energies using Eliashberg's technique |22| to obtain 
3xm(u) and 5xd{u). The zero frequency limit in Eq. 
(||) is performed analytically and the resulting equations 
are evaluated numerically. The sum over Q includes a 
summation over all Landau levels and over all dynamical 
fluctuation modes, in order to extend the theory beyond 
the region of validity for static, long-wavelength fluctua- 
tions, and beyond the lowest-Landau-level approach. 

The experimental zero-field transition temperature of 
92.5 K determines the temperature scale for the theoret- 
ical calculations. The mean-field transition temperature, 
T c (8.4 T) = 80.9 ±.3 K, which is determined by the diver- 
gence of the pair fluctuations, is obtained from our fit to 
spin susceptibility |10|| . We assumed H/2ttt^ — 0.02fcsT c 
and h/2TtT = 0.2fcsT c , and there is one fitting parame- 
ter for the overall scale of the fluctuation contributions 
to 1/Ti. Our theoretical calculation for the field depen- 
dence of the fluctuation correction is shown in Fig. ^ for 
d-wave pairing. The rate increases because of the sup- 
pression of the (negative) DOS contribution to the rate 
by the magnetic field. The results agree quantitatively 
with the experimental data at T = 95 K and provide 
strong evidence for d-wave pairing fluctuations. For s- 
wave pairing the calculated rate (not shown) decreases 
with increasing magnetic field because of the suppression 
of the (positive) MT term. 



3 




-0.25 — I 1 1 1 1 

5 10 15 20 25 30 

H - Field [T] 

FIG. 3. Fluctuation contribution R(l/(TiT)) = 
5(TiT) _1 / (TiT)n 1 of 63 Cu(2) spin-lattice relaxation rate 
as a function of magnetic field at 95 K. The dashed curves 
are d-wave calculations for temperatures ranging from 93 
to 102 K in increments of 1 K. The solid curve is calcu- 
lated for 95 K. The open circle is i?(l/(TiT))jvQfl at 95 
K [lisj. The open squares are from direct measurements 
of the 63 Cu(2) Ti at 3.5 T by Y.-Q. Song @, 5.9 T by 
Carretta et al. ||, and 7.4 T by Hammel et al. Q. 

Carretta et al. || reported experimental evidence for 
a positive contribution to the rate that was attributed 
to the MT process. These authors compared nuclear 
quadrupolar resonance (NQR) relaxation measurements 
and NMR relaxation at 5.9 T and found an NQR rate 
that is higher than the NMR rate in a range of ~ 10 K 
above T c , a result which is similar to our NQR measure- 
ment shown in Fig. |[ They interpret the decrease from 
the higher NQR to the lower NMR rate at 5.9 T in terms 
of s-wave pairing fluctuations, which implies a dominant 
MT term. However, our data in Fig. |^ shows that there 
is no significant MT contribution to the NMR rate at 
fields above 2.1 T. Our analysis of the field dependence 
of the data is in excellent agreement with the theory of d- 
wave pairing fluctuations, and disagrees with the theory 
based on s-wave fluctuations. Possible explanations for 
the apparent discrepancy between the NQR rate and the 
low-field NMR rate include an admixture of s-wave and 
d-wave fluctuations induced by orthorhombic anisotropy 
P), and the 2D to 3D crossover regime at low fields. 

In summary, we have determined the 63 Cu(2) spin- 
lattice relaxation rate as a function of magnetic field 
from 2.1 T to 27.3 T. We found that 1/TiT increases 



with increasing field in the temperature range T < 120 
K, which we can account for quantitatively with the the- 
ory of d-wave pairing fluctuations in 2D. Our results are 
consistent with d-wave pairing in YBCO, and inconsis- 
tent with dominant s-wave pairing. We found that the 
characteristic field scale for the suppression of the fluctu- 
ation corrections, 5(TiT) -1 , is ~ 10 T, which is an order 
of magnitude smaller than the expected field scale for a 
purely magnetic scenario for the pseudogap. 
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